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ABSTRACT 

We present the first unambiguous evidence of a broad (Gaussian width ~330 eV) 
component of the iron K a fluorescent emission line in the X-ray obscured Narrow 
Line Seyfert 1 Galaxy NGC 5506. This is the main results of a spectroscopic monitor- 
ing campaign on this source performed with the XMM-Newton observatory between 
February 2001 and January 2009. The broad line lacks extreme redwards skewness. If 
modelled with a relativistic component, the profile of the line is consistent with a flat 
emissivity radial dependence (a ~1.9). The disk inclination (~40°) is nominally larger 
then typically observed in unobscured AGN, in agreement with most measurements 
of broadened iron lines in Seyfert 2 galaxies. The quality of the data allows us to 
decompose the full iron emission line complex, and to study its long-term (timescales 
of weeks to years) variability pattern. The intensity of the neutral and narrow iron K Q 
core remains constant during the monitoring campaign. This indicates that the opti- 
cally thick gas responsible for the non-relativistic reprocessing of the primary AGN 
continuum in NGC 5506 is probably located in the torus rather than in the optical 
Broad Line Region. 

Key words: Galaxies: active - Galaxies:Seyferts - X-rays:galaxies - X- 
rays:individual:NGC 5506 



1 INTRODUCTION 

Broadening by relativistic effects of emission lines in an 
X-ray illuminated accretion disk offers the opportunity 
of probing the physical conditions and geometrical dis- 
tribution of matter in the vicinity of black holes, as 
well as General Relativity effects (Reynolds & Nowa k~2003[ 
Miller 2007). Broad and skewed profiles of the iron K a 
fluorescent emission line were discovered by ASCA in the 
mid-90s (|Tanaka et al. 19951 [Mushotzky et al. 1995[ ). High- 
throughput observations with the EPIC cameras on-board 
XMM-Newton have shown that at least 30 percent of 
nearby Seyfert galaxies exhibit relativistically blurred fea- 
tures consistent with being produced within 50 gravita- 
tional radii (r g = GM/c 2 , where M is the black hole mass) 
l|Nandra et al. 2007]) . However, doubts have been cast on 
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the robustness and uniqueness of these results. Ionised out- 
flows partially covering the primary high-energy emission in 
AGN could mimic relativistic effects ([Turner & Miller 2009 
I Reeves et a l. 2004). Controversy is fierce. In this context, it 
is crucial to confirm whether objects, whose data had poor 
signal-to-noise, do exhibit relativistically broadened lines as 
well as to expand the parameter space covered by these mea- 
surements. 

The study of X-ray obscured AGN is particu- 
larly promising. In the "Unified Model" framework 
(Antonucci & Miller 1985; Antonucci 1993), their putative 
accretion disks are expected to be observed at higher in- 
clination angles if the plane of the disk is aligned with 
that of the obscuring matter. For inclinations not too close 
to "edge-on", special relativistic effects, such as Doppler 
boosting, should be therefore stronger and easier to de- 
tect. Moreover, in type II objects the polar outflows should 
be out of the line of sight and this should break the de- 
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generacy between "warm absorbers" and relativistic lines 
i|Turner fc Miller 20091) . Finally relativistic lines allow the 
direct measurements of the inclination angle of the accre- 
tion disk thus offering a simple test of the unified scenarios. 

Observations of broad Fe line in X-ray obscured Seyferts 
are so far sparse. The Fe K a line in the Suzaku spectrum 
of MCG-5-23-16 implies a disk inclination, i, of 53° ±9° 
i|Reeves et al. 2007[PT . Nandra et al. (2007) add to this cen- 
sus NGC 526A (i=43°±fg°) and NGC 2992 (j=24°±7°), 
whereas NGC 2110 and Mkn 6 do not exhibit significant ev- 
idence for relativistic effects. NGC 1365 (i=24°±!°; Risaliti 
et al. 2009) and IRAS13197-1627 (i=27°±17°; Dadina & 
Cappi 2004, Miniutti et al. 2007b) were also recently dis- 
covered to host broad lines from nominally low inclination 
disks. The aforementioned uncertainties are purely statis- 
tics; systematic uncertainties related e.g. to different spec- 
tral deconvolutions are seldom discussed in the literature. 
They can be at least comparable to the statistical ones. 

In this context, the case of the nearby (z = 0.0061) 
X-ray obscured (N H ~ 3 x 10 22 cm" 2 ; Wang et al. 1999) 
Narrow Line Seyfert 1 Galaxy (NLSyl; Nagar et al. 2002) 
NGC 5506 is one of the most controversial. Bianchi et al. 
(2003) analysed all the X-ray observations performed with 
imaging and spectroscopic instruments on NGC 5506 at that 
time. They did not detect any spectral components directly 
related to the accretion disk. This led them to conclude that 
the disk is either fully ionised, or almost edge-on. Two early 
XMM-Newton observations of NGC 5506 (originally pre- 
sented by Matt et al. 2001) are also in the Nandra et al. 
(2007) sample. In their analysis, only one observation shows 
evidence of broadened emission. Nandra et al. (2007) inter- 
preted the lack of detection as simply due to low signal-to- 
noise ratio. 

In this paper we present spectroscopic results of a 
XMM-Newton monitoring campaign of NGC 5506, which 
brings the total integration time to ~197 ks (versus the 
23 ks of the observations published so far). The main re- 
sults of this paper is the first unambiguous detection of a 
broad iron K a line in NGC 5506. 



2 OBSERVATIONS AND DATA REDUCTIONS 

Tab. [1] lists the XMM-Newton observations of NGC 5506 
performed so far. We discuss in this paper data of 
the EPIC cameras only (MOS, Turner et al. 2001; pn, 
Struder et al. 2001). Data were reduced using the SASv9.0 
l)Gabriel et al. 2003P tasks e[mp]proc with standard set- 
tings. Periods of high-background were removed using count 
rate thresholds on the high-energy, single events, field-of- 
view light curves (Tab. [TJ which optimise the signal-to- 
noise of each individual observation. Source spectra were 
extracted from circular regions surrounding the source cen- 
troid (Tab. [TJ), whereas background spectra from circular 
regions on the same chip free from contaminating serendip- 
itous sources (and at the same height in detector coordi- 
nates for the pn to ensure that the same Charge Transfer 
Inefficiency correction applies). Response files were gener- 
ated with the SAS tasks arf gen and rmfgen. Spectra were 

1 see, however, Nandra et al. (2007), whose spectral deconvolu- 
tion of the XMM-Newton/EPIC spectra requires 1^20° 



Table 1. Log of XMM-Newton NGC 5506 observations 



XMM-N ID 


Start 


P a 


T C R b 


- 1 - exp 




Time 


(") 


(s- 1 ) 


(ks) 


0013140101 (l) d 


02-02-2001 


100/40 


0.2/0.5 


13.3 


0013140201 (2) d 


09-01-2002 


120/44 


0.5/1.0 


9.9 


0201830201 (3) d 


11-07-2004 


42/45 


0.5/0.5 


14.8 


0201830301 (4) d 


14-07-2004 


40/42 


0.5/0.35 


13.9 


0201830401 (5) d 


22-07-2004 


40/42 


0.35/0.35 


13.8 


0201830501 (6) d 


07-08-2004 


40/42 


0.35/0.35 


13.9 


0554170101 (7) d 


27-07-2008 


50/50 


0.35/0.35 


55.8 


0554170201 (8) d 


02-01-2009 


40/40 


0.5/0.35 


62.0 



"size of the source spectrum extraction region in the MOS/pn; 
6 count rate threshold applied to select intervals of low particle 
background in the MOS/pn; c pn exposure time; d label used to 
identify the observation in this paper: XMMi 



rebinned to over-sample the intrinsic energy resolution by 
a factor not larger than 3, and to ensure that each spec- 
tral channel has a number of background-subtracted counts 
larger than 50. Spectral fits were performed in the 2.2- 
10 keV energy band, to avoid large gradients of the instru- 
ment transfer function due to the instrumental edges (1.8- 
2.2 keV) and to the soft X-ray photoelectric cut-off. Below 
the energy of the latter the X-ray spectrum is dominated 
by reprocessing and scattering of the nuclear radiation by 
photoionised gas on scales as large as a few hundreds par- 
sees (|Bianchi et al. 2003[) . High-resolution spectroscopy of 
this component will be discussed elsewhere (Labiano et al., 
in preparation). 

In this paper: energies are quoted in the source rest 
frame; errors are quoted at the 90% confidence level for 
one interesting parameter; the following cosmological pa- 
rameters have been adopted in the calculation of the lu- 
minosities: #o=70 km s" 1 Mpc" 1 , A o =0.73, fi m =0.27 
i|Bennett et al. 2003|) . 



3 SPECTRAL ANALYSIS 

3.1 Phenomenological fits on individual 
observations 

As first step, we fit the pn spectrum of each individual 
XMM-Newton observation with the same phenomenologi- 
cal model. Following the approach as in Matt et al. (2001), 
the model is constituted by a photoelectrically absorbed 
power-law, a photoelectric absorption edge and three nar- 
row iron emission lines, modelling the K Q fluorescence of 
Fel and of recombination lines of Fexxv and Fexxvi. The 
Fel and Fexxvi are actually each close doublets, which can- 
not be resolved by the EPIC cameras. Their centroid en- 
ergy has been fixed to the weighed mean of the correspond- 
ing doublet components: 6.399 keV and 6.966 keV, respec- 
tively. The structure of the Fexxv line is instead more com- 
plex, with the resonance (6.7002 keV), two intercombination 
(6.6821 keV, and 6.6673 keV) and the forbidden components 
(6.6364 keV) all potentially contributing. In the fit of each 
observation we left originally the Fexxv centroid energy free 
to vary. We obtained lower limits consistent with the en- 
ergy of the resonance line (and inconsistent with the other 
components) in all observations, except in XMM8, where 
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Figure 1. 5.5-7.5 keV residuals in units of standard deviations 
when the best-fit phenomenological model is applied to the 2.2- 
12 keV EPIC-pn spectra of the 2004-2009 XMM-Newton obser- 
vations of NGC 5506. 



E FeXXV = 6 654± ; 020 ke y = (3 7 ± LQ ) x 1Q -5 pho _ 

tons cm -2 s -1 ]. In all the subsequent models we fixed the 
centroid energy of the Fexxv line to the value of the reso- 
nance component, except for XMM8 for which we used the 
value of the forbidden component. Finally, the energy of the 
edge was constrained to be <;7.11 keV, corresponding to 
photoionisation of iron. 

The Fel line 1-a width (only "width" hereafter) is 
46±J4 eV. This measurement is consistent with the width 
of the Manganese K a and Kp pn calibration emission lines 
measured in observations close to the 2009 NGC 5506 ob- 
servations: a ca i = 42 ± 8 eV, and a ca i = 47 ± 8 eV in 
Obs.#0410780601 (August 9, 2008) and Obs.#0412580401 
(January 4, 2009), respectively. The measurement is thus 
consistent with the instrumental resolution. We have con- 
sequently used this value for all the nominally unresolved 
astrophysical emission lines in the models discussed in this 
paper. 

In Fig.[T]we show the residuals against the best-fit phe- 
nomenological model in the 5.5-7.5 keV energy band (the 
residuals are flat outside this interval) . In the longest obser- 
vations significant wave-like residuals are present around the 
nominal energy of the Fel fluorescent K„ line. Hints of fea- 
tures with a similar shape are visible in other observations 
with a lower statistics. These features can be interpreted as 
being due to a broad (a ~300 eV) excess emission feature, 
which the narrow Fel Gaussian profile is unable to account 
for. 

In Fig. [T] we do not show the residuals for the first two 
observations (XMM1 and XMM2). The event pattern dis- 
tribution measured by the SAS task epatplot suggests that 
they might be marginally affected by pile-up due to the dif- 
ferent instrument mode (Large Window instead of Small 
Window). This may also explain why XMM2 - which corre- 
sponds to the highest X-ray flux measured during our cam- 
paign - exhibits the flattest X-ray spectral index. This is 
illustrated in Fig. [2] where the best-fit parameters of the 
phenomenological model in each individual fit as a function 
of time and of the 7-12 keV absorption-corrected flux are 
shown, once the excess emission around 6 keV is fit with a 
further Gaussian profile whose width was left free to vary. 



The intensity of the narrow Fel line is constant. Its mean 
normalised dynamical range is ^ 27%. Evidence for variabil- 
ity of the recombination He- and H-like line intensity is also 
marginal. Beside the possible change of ionisation state of 
the Fexxv in XMM8, it is limited to a nominal difference 
by a factor of at least five between two measurements of the 
Fexxvi line during the 2004 monitoring campaign. Given 
the uncertainties in the spectral deconvolution at ~7 keV, 
we consider this evidence for variability as only tentative, to 
be confirmed by high-resolution measurements. 

We have applied the same phenomenological model 
also to the combined spectra extracted from the MOS 
cameras. The parameters of the continuum are different. 
The ratio between the 2-10 keV MOS and pn fluxes 
(weighted mean, (F M os / F pn )=5.5%; o FMOS / Fpn =2.7%) is 
consistent with known systematic calibration uncertainties 
in the absolute flux determination (Stuhlingcr e t al. 2008] ). 
However, the average difference between the spectral 
index ((Ar)=0.23; o"Ar=0.07) and the column density 
((AiVH)=4.9xl0 21 cm" 2 ; a A N H =2. Ox 10 21 cm" 2 ) is larger 
then expected. The difference is larger for higher fluxes. The 
continuum spectral parameters measured by the EPIC cam- 
eras are instead well consistent for observations whose 2- 
10 keV flux is lower than ~8xl0~ n erg cm" 2 s" 1 . Visual 
inspection of the MOS pattern fraction distribution unveiled 
an excess (deficit) of double (single) events by 15-30% above 
4 keV in higher flux observations. The above pieces of ev- 
idence indicate that many of the NGC 5506 MOS spectra 
could be still affected by pile-up. Consequently, we will not 
discuss the MOS spectra in the rest of the paper. However, 
we will show in Sect. 3.2 that the main result of this paper 
- the ultimate discovery of a broad component of the iron 
K„ line - is measured in both the pn and the MOS cameras 
with comparable properties. 



3.2 Physically-motivated combined fits of all 
observation spectra 

3.2.1 Description of the models 

Motivated by the results presented in the previous Sec- 
tion, we enriched the phenomenological model by adding 
physically-motivated components to account for features de- 
tected in the spectral residuals. We have investigated three 
main scenarios: 

• Scenario #1: relativistically broadened K a iron line 
and non-relativistic Compton reflection. The former com- 
ponent has been modelled through the code kyrline 
UDovciak et al. 2004[l . which calculates the profile of an 
emission line produced in an axisymmetric accretion disk 
around a black hole. The Compton-reflection component 
(model pexrav; Magdziarz & Zdziarski 1995) is primarily 
associated in this scenario to the narrow core of the iron K a 
line, although it may include a contribution by the accre- 
tion disk as well, which the moderate energy resolution of 
the EPIC cameras is unable to disentangle. 

• Scenario #2: non-relativistically broadened K a iron line 
and Compton reflection. The former component has been 
modelled with a simple Gaussian profile 

• Scenario #3: relativistically broadened emission line 
and continuum reflection dominated by a relativistically 
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Figure 2. Best-fit parameters of the phenomenological fits of individual NGC 5506 XMM-Newton observations as a function of time 
(left panel) and of the 7-12 keV absorption-corrected flux (right panel): from top to bottom: photon index, column density, intensity 
of the Fei, Fexxv, FeXXVI narrow emission lines and of the broad Gaussian. Continuum measurement data points corresponding to pn 
observations in Large Window mode are surrounded by a large circle. 



blurred component. The latter component has been mod- 
elled by convolving pexrav with the relativistic kernel 
kyconv ( Dovcia k et al. 2004fl and assuming a constant ratio 
between the normalisation of the primary continuum and 
that of the Compton reflection (R, see Tab.(2])- 

kyrline belongs to a suite of models based on a com- 
mon ray-tracing subroutine aiming at describing the X-ray 
emission of black-hole accretion disks in the strong gravity 
regime ( Dovciak et al. 2004 ). The line profile depends on the 
following parameters: a) the black hole spin a (comprised 
between and 0.9982 in dimensionless relativistic units); b) 
the inner (ri n ) and outer (r out ) radius of an annular region 
on the disk where the line photons are emitted; c) the index 
q of the radial dependence of the emissivity per unit area, 
ft, in the local frame comoving with the disk: «(r) oc r~ a ; d) 
the "disk inclination" i, i.e. the angle between the normal 
to the disk plane and the line-of-sight; e) the rest frame en- 
ergy (assumed monochromatic) of the photons emitted by 
the disk, E c . Although the line normalisation in kyrline 
is expressed in units of photons cm -2 s _1 integrated over 
the whole profile, we will also use the line Equivalent Width 
(EW), i.e. the line intensity normalised to the underlying 
continuum at 6.4 keV, to ease comparison against theoreti- 
cal predictions IMatt et al. 1992]) . The convolution function 
kyconv uses the same relativistic kernel as kyrline. 

The continuum reflection components replace the 
photo- absorption edge in the phenomenological model. We 
always linked together the spectral index of the Compton- 
reflected continuum and that of the primary emission. More- 
over, we have added two narrow-band spectral compo- 
nents to the non-relativistic emission line complex, mod- 
elling the: a) iron K a fluorescent iron line Compton-shoulder 
( |Sunyaev fc Churazov 1996[ IMatt 2002]) through a rectan- 



gular box function comprised between 6.24 and 6.4 keV, 
and whose intensity was constrained to be not larger then 
20% the intensity of the corresponding narrow K a line (the 
fits determine only an upper limit on this parameter as 
large as the maximum allowed value); b) the Kp line asso- 
ciated with the iron K a narrow component, modelled with 
a Gaussian profile with energy fixed at 7.058 keV and in- 
tensity constrained not to exceed 16% of the K a intensity 
HMolendi et al. 2003[1 . 

In principle, only model parameters describing physical 
quantities not expected to change among the observations 
(or demonstrated to remain constant in the phenomenolog- 
ical analysis) have been constrained to have the same value 
for all fitted spectra. They are: the intensities of the narrow 
Fei, Fexxv and Fexxvi lines; the metallicity (assumed so- 
lar; Anders & Grevesse 1989); the high-energy cut-off of the 
intrinsic primary continuum (held fixed to 130 keV in accor- 
dance to the results of the BeppoSAX observation; Bianchi 
et al. 2003); the inclination of the disk and of the distant 
non-relativistic reflector (the latter held fixed to 45°); and 
the black hole spin. In Scenario#l and #3 varying the black 
hole spin from Schwarzschild to maximally rotating yields a 
variation of the \ 2 lower than 0.5. We have therefore fixed 
this parameter to the value corresponding to the latter so- 
lution (a — 0.998), favoured at face value. Once these as- 
sumptions were made, a first series of fits showed that many 
of the remaining free parameters were consistent with being 
constant across all the fitted observations. We have there- 
fore forced them to assume the same value in all the fits. 
In Scenarios #1 and #3 they are: a) the inner and outer 
radius of the X-ray emitting disk. Initially we have fixed the 
latter to 400r g , because this parameter is generally uncon- 
strained when the line profile is strongly relativistic. We will 
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Figure 3. Residuals in units of standard deviation when the Sce- 
nario#3 best-fit model is applied to the 2004-2009 EPIC-pn spec- 
tra. From top to bottom: XMM3, XMM4, XMM5, XMM6, XMM7, 
XMM8. Y-axis range and colours are the same as in Fig. [T] 



review this assumption in Sect. 3.2.2. The former has been 
constrained throughout the paper to coincide with the in- 
nermost stable circular orbit of an accretion disk co-rotating 
with a black hole of spin a; b) the index of the power-law 
disk emissivity radial dependence, a; c) the energy of the 
relativistically broadened line in the source rest frame, E c . 
In Scenario #2 they are all the parameters describing the 
broad Gaussian profile. 



3.2.2 Results 

The Scenarios described above yield comparable fit quali- 
ties: x 2 /^= 12 65. 3/1161, 1280.9/1161, and 1258.1/1161 for 
Scenario #1, #2, and #3, respectively. These x 2 are signifi- 
cantly better than those yielded by the same models without 
the broad line: 1479.5/1169, 1516.2/1169, and 1581.6/1166, 
respectively. No further structure is visible in the residuals 
(Fig. [3]). We therefore concluded that no additional astro- 
physical component is required. In Tab. [2] we summarise the 
best-fit results for the non-relativistic components in Sce- 
nario^. The values for the other scenarios are consistent 
with them within the statistical uncertainties. In Fig. [4] we 
show the spectra of observations XMM3 to XMM8 together 
with the components of their best-fit models. 

The properties of the relativistic line can be summarised 
as follows (see Tab. [3]): a) the rest frame centroid en- 
ergy is consistent with neutral or moderately ionised iron 
(E c ~ 6.48 keV). Fixing the source frame centroid energy to 
the value expected for neutral iron in, e.g., Scenario#3 yields 
X 2 =1276.0/1160 dof; b) the disk is slightly more inclined 
(l ~ 40°) than typically measured in unobscured Seyfert 
galaxies (|Nandra et al. 19971 INandra et al. 2007|) . and; c) 
the line profile is only mildly relativistic as parametrised 
by the moderately flat (a ~1.9) radial emissivity profile. 
The line emitting region extends close to the innermost sta- 
ble orbit. The innermost radius is constrained between 3 
and 25 r g . In Fig. [5] we show the profile of the relativisti- 
cally broadened iron line, obtained by stacking together the 
source frame residuals of each individual observation against 
the best-fit model in Scenario#3 once the relativistically 



Table 2. Best-fit results for the parameters of the non-relativistic 
continuum components in Scenario #3 (except parameter EWp e j 
which refers to Scenario#l). 



Spectrum dependent parameters 
Obs.# 



22 



2-10 keV Flux a 



XMM3 
XMM4 
XMM5 
XMM6 
XMM7 
XMM8 



(10 
2.86± 



cm ") 
0.11 
0.14 
2.94±0;l| 

2.65+0;!! 

2.99±0;?» 

3.08±°;0| 



1.778± 



_0.014 
-0.025 
1 777+0-D13 

±- ' ' ' ^0.033 

i 771+0.014 

1.844±°;° 1 j 
1.864 ± 0.009 

J). 010 
-0.006 



(io- 



82.5±°/f> 
79.4±i : f 

120.0±?J 
132. 5±H 



1.862±K 



118. 7± 



0.7 



Spectrum independent parameters 
(5.9 ±0.6) X 10~ 5 °-" -™" 2 



jb 

1 FeXXVr 



erg cm 

1.0 X 10~ b erg cm~ 2 

r i —9 



I0 FeXXVf ^ 14 X 10 " er S Cm 
4eXXV/=( 2 -°±°-7)xl0 

R s; 0.1 c 

EW FeI > 1.0 keV d 



s 

erg cm" 



"corrected for absorption. XMM1: 70.8 ± 1.0; XMM2: I2\.l±\\ 
in the same units; ^intensity of the narrow emission lines; c ratio 
between the primary continuum and the relativistically blurred 
reflection normalisation; d Equivalent Width of the narrow com- 
ponent of the Fel fluorescent line against its own reflection con- 
tinuum. 

Table 3. Best-fit values of the K a iron relativistically broadened 
profile in Scenarios #1 and #3 for r oll t=400r 9 . The rightmost 
column contains the statistical uncertainties to be added if r ou t 
is allowed to vary. 







Scenario #1 


Scenario #3 


Additional error 
if r ou t is left free 




(keV) 


6.48 ±0.02 


6.48 ±0.02 


^0.02 
^0.00 


a 




1.9 ±0.3 


1.9 ±0.3 


_i_0.2 
^O.l 


i° 




42±<! 


37 ±5 


±r 


' in 


ir g ) 


15±? 2 


15±J1 
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q 
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c 

o 
o 
15 
IS 
Q 

1.00 



0.95 



NGC5506 



^ ******* 



6 7 8 

Energy (keV) 



Figure 5. Stacked profile of the relativistically broadened iron 
line. The dashed line indicates the observed frame energy of the 
neutral iron K a fluorescent line. 
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Figure 4. EPIC-pn NGC 5506 spectra in observations XMM3 to XMM8 (crosses). The lines indicate different model components: 
dot-dashed: continuum; long dot-dashed: narrow lines; solid: broad line. 



broadened emission line had is removed from the model. The 
profile is not significantly skewed. This is consistent with 
the similar fit quality that relativistic and non-relativistic 
scenarios yield. A small "dip" in the profile coincide exactly 
with the rest frame energy of the K a neutral iron fluorescent 
line. It could be due to an incorrect disentangling between 
the narrow and the broad profile, whereby the fit preferen- 
tially attributes photons in this energy range to the former. 



We estimate that the intensity of the narrow component 
should be lowered by at least 50% for the dip to disappear. 
The effect of this potential systematic uncertainty on the de- 
termination of the relativistic line profiles a and Vi n is negli- 
gible, because they are mostly sensitive to the red wing. On 
the other hand, centroid energy and inclination angle could 
be affected. Assuming If b i = 3 x 10 -5 photons cm -2 s _1 , 
they change by AE C ~ 30 eV and At ~ 4°, respectively. 
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Figure 6. Iso-x 2 contour plots for the intrinsic width of a Gaus- 
sian profile fitting the broad component of the iron K a line in the 
NGC 5506 merged spectra. The contours correspond to the 68%, 
90% and 99% confidence level for two interesting parameters, re- 
spectively. Solid lines: pn versus MOS1; dashed lines: pn versus 
MOS2. 

Nonetheless, the \ 2 significantly worsens (1361.9/1160 dof 
in Scenario#3). This hypothesis is then inconsistent with 
the data. 

We come back now to the original restriction imposed 
to the value of r out . The flat value of a implies that a not 
negligible contribution to the broad line profile comes also 
from large disk radii. We have therefore repeated the fits in 
Scenario#3 after allowing r out to be a free parameter. At 
the 90% confidence level for one interesting parameter r ou t 
is constrained to be comprised between 250 and 900r 9 . The 
best-fit value is still around 400r 9 , hence the best-fit values 
do not change. However the broad line parameters errors are 
augmented by a small amount (Tab. [3]). 

In Scenario #2, the centroid energy of the Gaussian 
profile is consistent with fluorescence from neutral or mod- 
erately ionised iron: E c — 6.45±o;ol keV. The line width is 
o h = 330 ± 40 eV. 

We have checked the robustness of the broad iron 
line profile against possible data reduction and/or undi- 
agnosed calibration problems. A fit with Scenario#2 on a 
merged pn spectrum with single and double events yields 
the following best-fit broad Gaussian profile parameters: 
E c = 6.4T±g;8t keV, a b = 360±|g eV, and I b = (1.47 ± 
0.14) x 10 -4 ph s _1 cm* 2 . They agree within the statisti- 
cal uncertainties with the best-fit parameters obtained on 
a merged pn spectrum extracted with single events only: 
E c = 6.41 ± 0.04 keV, a b = 360±|g eV, and I b = (1.55 ± 
0.18) xlO -4 ph s _1 cm . More importantly, the width mea- 
sured by the pn spectrum agrees with that simultaneously 
measured by the MOS cameras (Fig. [6]). 

3.2.3 Variability 

The EW of the broad component of the iron K a line is con- 
sistent with being constant among the XMM-Newton obser- 
vations discussed in this paper (Tab. [4]). On the other hand, 
its intensity is not constant. A fit on the broad line intensity 
versus time relation with a constant yields a x 2 /" / =H-4/5 
(versus 2.1/5 for the EW). This indicates that in the first 
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Table 4. Measurements of the broad iron K a line EW (in eV). 



Obs.# 


Scenario #1 


Scenario #2 


Scenario #3 


XMM3 
XMM4 
XMM5 
XMM6 
XMM7 
XMM8 


130±2g 
180±f° 
140 ± 30 
122±i 
152±" 

149±i| 


130 ± 30 
160 ± 30 
130 ± 30 

ioo±lg 

137±}3 
143±» 


120±18 
150±3° 
120 ± 30 
110 ±20 
132±J3 
131 ± 15 


XMM78i 
XMM78 2 
XMM78 3 
XMM78 4 


180±§g 

i±u=c 20 
120 ± 20 


140±1° 

90 ±20 
100 ± 20 


160 ± 20 

100±lo 
101±« 

109±« 



approximation the broad line intensity follows the variation 
of the underlying continuum. 

The absorption-corrected 2-10 keV flux measured dur- 
ing the XMM-Newton observations covers a dynamical range 
~1.9. This is less than a factor of two smaller than mea- 
sured during the 6-years long intense X-ray monitoring with 
the RXTE/PCA in the same energy band (~3.5; Uttley & 
McHardy 2005). The longest time difference between two ob- 
servations, whose continuum flux differs by more then 27% 
(the upper limit on the percentage fractional variability of 
the intensity of the iron K Q fluorescent line narrow compo- 
nent) is ~7.92 years (XMM8 versus XMM1). 

The photon index in observations XMM6 to XMM8 
(2-10 keV absorption corrected flux, F2-10 5s 1.19 x 
10 -10 erg s _1 cm -2 ) is systematically steeper than in the 
others (by Ar ~0. 08-0. 10). There is an intrinsic correla- 
tion between the power-law spectral index and the column 
density of a photoelectric absorber covering it as measured 
by the spectral fitting. Steep spectra invariably correspond 
to higher obscuration. In order to disentangle the intrinsic 
driver of the continuum spectral variability we run two Sce- 
nario's models, where the power-law index or the column 
density had been constrained to assume alternatively the 
same value among all spectra. The quality of the fit - albeit 
worse than when both parameters are allowed to vary - is 
better in the case when the power-law is allowed to vary: 
X 2 ' /v — 1276.4 against 1300.6/1164, suggesting that varia- 
tions of the intrinsic continuum are the main driver of this 
long-term continuum spectral variability. 

The X-ray light curve of NGC 5506 shows 
also rapid variability, by a factor ~1.7 over 250 s 
(ALx==;l.lxl0 42 erg s _1 ; Dewangan & Griffiths 2005). 
The variability is achromatic (Fig. [7] McHardy & Czerny 
1987). In order to investigate the possible dependence 
of the broad line parameters upon the overall source 
flux on such small timescales, we have combined the 
event lists of the 2009 observations, and extracted four 
intensity-resolved spectra, corresponding to the following 
background-subtracted count rate ranges: <2.9, 2.9-3.2, 
3.2-3.5, and >3.5 s _1 . These ranges were chosen in such 
a way that the intensity-resolved spectra have approxi- 
mately the same net background-subtracted counts. These 
spectra are labelled as XMM78j in Tab. |4] with j=l to 
4 (in increasing flux order). The corresponding 2-10 keV 
absorption-corrected fluxes are: i09.3±i:!, i24.3±i:l, 
134.6±i:l, 145.0±i;?xl0~ 12 erg cm -2 s" 1 , respectively. No 
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Figure 7. 3-5 keV versus 1.5-3 keV hardness ratio as a function 
of the 1.5-5 keV count rate in all the XMM-Newton observations 
of Tab. [T] 
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Figure 8. Residuals (in units of data/model ratio) when the best- 
fit continuum in Scenario#3 is applied to intensity-resolved spec- 
tra in the nominal energy range excised of the 5-7 keV interval. 
The contribution of the non-relativistic emission line components 
was removed as well by subtracting their time-averaged intensity 
from the residuals. Filled circles, XMM78i; crosses, XMM782; 
empty circles, XMM783; stars, XMM784. 

correlation exists between the EW of the broad K Q iron 
line component and the flux (Tab. 0J, save for a larger EW 
in the spectrum corresponding to the lowest flux state. The 
shape of the broad line component remains also remarkably 
constant in different flux states (Fig. [8} 



4 DISCUSSION 

4.1 Continuum variability 

NGC 5506 is one of the most rapidly variable AGN of the 
X-ray sky. It exhibits luminosity changes of ~10 42 erg s _1 
within a few minutes (Dewangan & Griffiths 2005). The 
variability is, however, remarkably achromatic on these time 
scales (Fig. see also McHardy & Czerny 1987). On longer 
timescales, a small but systematic difference in the intrinsic 
spectral shape (Ar ~ 0.1) is found between "low flux" and 
"high flux" observations (absorption-corrected flux ratio be- 



tween the latter and the former of about 1.5). An even larger 
dynamical range in both luminosity and spectral shape was 
measured during a hard X-rays (E > 20 keV) monitoring 
campaign with HEXIT (Manchanda 2006). The variability 
pattern measured by XMM-Newton is inconsistent with a 
pair-dominated corona (Haardt ct al. 1997), where an in- 
crease of the spectral index by the amount observed by 
XMM-Newton requires at least an order-of-magnitude in- 
crease of the 2-10 keV luminosity. This does not automati- 
cally imply that pairs are not an important source of scatter- 
ing in the Comptonizing material; the alternative possibility 
is that the increase of the corona optical depth is due to a 
variation of the linear dimension of the corona. However, the 
shape of the broad iron K a iron line does not change signif- 
icantly on the shortest possible variability timescales which 
XMM-Newton can probe (Fig. [8]), suggesting that the illu- 
mination of the disk does not significantly change. 

4.2 The origin of the narrow iron K Q line 

One of the main results presented in this paper is that the 
intensity of the narrow component of the K Q in NGC 5506 
remains constant during the whole eight years covered by 
the XMM-Newton campaign. This extends by a factor of 
3 the longest baseline published so far - a ~1000 days 
RXTE monitoring campaign at the end of last century 
l|Lamer et al. 2000|) . The RXTE/PCA spectroscopic mea- 
surements were unable to disentangle the two components 
(narrow and broad) of the iron K a line. Nonetheless Lamer 
et al. (2000) claimed that the strength of the contin- 
uum reflection component and the total EW of the iron 
lines were stronger during phases of low flux. The origin 
of the iron K a emission line has challenged the ingenu- 
ity of observational and theoretical astronomers since its 
early discovery ( |Mushotzky 1982| ITurner fc Pounds 1989)) . 
Evidence that this spectral feature is almost invariably 
accompanied by a Compton reflection continuum compo- 
nent ([Pounds et al. 19901 INandra fc Pounds 19"94)) in un- 
obscured AGN suggested an origin in optically thick ma- 
terial out of the line of sight. The origin of very broad 
and skewed profiles in the innermost region of the accre- 
tion disk is a matter of heated debate, that we briefly 
summarised in Sect. 1 (we will further discuss this point 
in Sect 4.3). As far as narrower (typical Full Width Half 
Maximum, FWHM, ,$3000 km s _1 ) and symmetric lines 
are concerned, possible physical locations for the line- 
emitting matter are the external regions of the accretion 
disk (Fabian ct al. 1989), the Broad Line Regions (BLRs) 
( |Yaqoob et al. 2 001 ; Bianchi ct al. 2008)) or the molecular 
torus l|Ghisellini et al. 19941 IKrolik et al. 1994)) . These sys- 
tems have sizes from a fraction to several tens of parsecs. 
The lack of correlation between the FWHM of the iron K a 
and that of the H^ lines (or the black hole mass) disfavours 
a general origin in the BLRs (jNandra 2006 [) . 

We can use the lack of variability of the narrow com- 
ponent of the iron K a fluorescent line to estimate the di- 
mension of its production region. The line produced by a 
reflector at distance d from an illuminating source tracks the 
variability pattern of the primary continuum if d ^ t var /c, 
where t vaT is the continuum variability timescale and c is the 
speed of light. We estimate t var as the longest time inter- 
val between XMM-Newton observations, whose 2-10 keV 
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absorption-corrected fluxes differ by more than 27%, i.e. 
than the upper limit on the percentage fractional variability 
of the iron line intensity. t va r hence corresponds to the span 
between observations XMM1 and XMM8: 7.92 years. This 
estimate is qualitatively in agreement with the determina- 
tion of the Power Spectrum Density by the RXTE/PCA 
monitoring campaign (Uttley & McHardy 2005), which ex- 
hibits a slope of ~-l below v ~ 4 x 10 -5 Hz and down to 
~6xl0~ 9 Hz. The XMM-Newton results constrain the lo- 
cation of the optically-thick matter responsible for the bulk 
of the iron K a narrow component at a distance <;2.5 pc. 
This is consistent with the best estimate of the inner side 
of the torus (< 5 pc; Prieto & Meisenheimer 2004). Unfor- 
tunately no direct estimate of the Broad Line Region size 
is available for NGC 5506. The 2-10 keV X-ray luminosity 
from the data presented in this paper is in the range ~6- 
llxlO 42 erg s _1 , which corresponds to a typical size of the 
BLR ^10 days (Kaspi ct al. 2005). This makes the torus the 
most likely origin of the narrow iron K a line. 

We report a possible change in the intensity of the 
Fexxvi during two of the observations in 2004 performed 
a few weeks apart, as well as a possible ionisation change 
of the Fexxv (from a resonance- to a forbidden-dominated 
state) during XMM8. Given, however, the spectral complex- 
ity around 7 keV, a confirmation of these findings through 
instrumentation with a much better energy resolution is 
needed. This is shown, for instance, by the different results 
on the forbidden component of the He-like during XMM8 
when one compares the phenomenological (Sect. 3.2) and 
the relativistic model (Tab. [2]). 

4.3 The origin of the broad iron Ka line 

This paper reports the discovery of a broad (a ~330 eV) 
component of the iron K a fluorescent emission line in 
NGC 5506. The line profile lacks a prominent red wing, in 
contrast to what typically observed in unobscured Seyfert 
galaxies (Miller 2007). Although fits with a relativistically 
broadened profile yield marginally better statistical quality 
than symmetric Gaussian broadening, it is worth posing the 
question of whether this difference may be due to a different 
physical mechanism responsible for the observed broaden- 
ing. 

The symmetric and highly localised profile rules out 
an explanation in terms of blending of transitions corre- 
sponding to different iron ionisation states, enhanced con- 
tinuum reflection (Guainazzi 2002; Risaliti 2002), scattering 
from a relativistic wind (Titarchuck ct al. 2009) or incor- 
rect fitting of the opacity affecting the underlying contin- 
uum (Turner & Miller 2009). Comptonization of line pho- 
tons was suggested to explain the first observations of broad 
lines by Ginga (Czcrny ct al. 1991). The observed moderate 
width requires the line photons to cross at most 3 Comp- 
ton depths of cool (kT c < 0.7 keV; a = y/2kT/mc 2 ) gas 
( [Pozdnyakov et al. 1979) |Sunyaev fc Titarchuk 1980") . The 
main shortcoming of this scenario, however, is the non- 
linear dependence of the line EW on the continuum flux 
on timescales as short as a few hundreds seconds, which is 
at odds with the line-emitting region being screened by op- 
tically thick matter. 

The explanation in terms of a relativistically broad- 
ened and skewed profile is consistent with the data. For- 



mal fits require a moderate inclination (i ~40°) accre- 
tion disk. NGC 5506 joins the so far restricted club of 
obscured AGN, which require relativistic effects from a 
moderately inclined (i ~ 40°) accretion disk, together 
with MCG-5-23-16 (IBalestra et al. 2004P and NGC 526A 
l|Nandra et al. 2 007 ). All the above sources exhibit a col- 
umn density obscuring the active nucleus in the range 
10 22-23 cm -2 . Although it is tempting to associate these 
sources with AGN seen through the rim of the obscur- 
ing torus (see, for instance, the discussion in Matt et al. 
2003), it should be borne in mind that several unobscured 
AGN in the Nandra et al. (2007) also formally require 
large inclination angles (e.g. Akn 120, Mkn 590, NGC 7213, 
NGC7469). On the other hand, the recently discovered 
broad iron lines in IR AS13197-1627 flDadina fc Cappi 2004 
IMiniutti et al. 2007aP and NGC 1365 (|Risaliti et al. 200*9 
require moderate disk inclinations of t ~ 27° and 24°, 
respectively. A comparison between unbiased samples is 
needed before any inference on the geometrical configura- 
tion of the torus and the disk in X-ray obscured AGN can be 
drawn from X-ray spectroscopy of the relativistically broad- 
ened K a line. Moreover, the disk inclination in NGC 5506 is 
significantly smaller than the inclination of the galaxy disk 
i|Imanishi 2000p . 

The broad line is not extremely relativistic. It lacks 
the prominent red wing extending down to ~4 keV 
observed in, e.g., MCG-6-30-15 ( |Fabian fc Vaughan 2003] 
IMiniutti et al. 2007ap , which requires very steep radial emis- 
sivity profiles (I Wilms et al. 200"T|> . The inner radius of the 
line emitting region (constrained by the data to be com- 
prised between 3 and 25 r 9 ) and the flat emissivity pro- 
file suggest that a large area on the disk contributes to the 
average profile. The line profile is remarkably constant in 
intensity-resolved spectra, which sample variations of the 
intrinsic flux on time scales as low as a few hundred seconds 
(cf. Fig. [5]) (against at odd with MCG-6-30-15; Iwasawa et al. 
1999) . The fits do not require statistically significant changes 
in the structure of the innermost accretion flow. The param- 
eters a and r; n describing the profile of the relativistically 
broadened iron line are consistent with being constant over 
5 epochs of observations spanning 9 years. 

Little is known of the black hole mass in NGC 5506. 
Both the central stellar velocity dispersion (~180 
km s" 1 ) l|Qliva et al. 19991 IGebhardt et al. 20001 
IFerrarese fc Merrit 20001 |Papadakis 2004) and the width 
of the [OIII] line l|Boroson 2003|> suggest a black hole 
mass ~10 s Mq. For this mass the shortest measured X-ray 
variability timescale implies that the bulk of high-energy 
continuum production and reprocessing occur within a few 
Schwarzschild radii from the black hole. This is in disagree- 
ment with the lack of extreme relativistic line broadening, 
as well as more generally with the idea that NLSyl Galax- 
ies are characterised by small black hole masses and high 
accretion rates (^2% in this scenario) (Koni 6ssT200"7)) . 
However, X-ray-based estimates suggest that the black hole 
mass in NGC 5506 could be significantly lower. McHardy 
et al. (2006) observe that lowering the black hole mass by 
a factor ~5 would bring NGC 5506 to the locus in the 
PDS break frequency, black hole mass and accretion rate 
plane shared by AGN and Galactic Black Hole systems. 
An even more radical revision (Mbh — 5 x 10 6 M©) has 
been proposed by Nikolajuk et al. (2009) on the basis of 



10 M. Guainazzi et al. 



the correlation with the X-ray 2-10 keV light curve excess 
variance. A black hole mass of this order (~2xl0 6 M Q ) 
had been proposed by Hayashida et al. (1998) using 
similar arguments. The light crossing time of the shortest 
variability timescale if one assumes the Nikolajuk et al. 
(2009) estimate corresponds to about 20 Schwarzschild 
radii, in better agreement with broad iron line spectroscopy 
results. In this case the accretion rate is ~40% of the 
Eddington limit, if one estimates the bolometric luminosity 
by applying a luminosity-dependent bolometric correction 
to the 2-10 keV luminosity (|Marconi et al. 2004[) . The 
size of the BLR, tblr ~ GM/v 2 , is ~8 days, if one 
estimates v from the F WHM of the Pa/3 broad component 
( |Nagar et al. 2002[ ) . 

The size of the AGN sample with robust detections of 
relativistic lines is still far too small for correlations with 
other observables to yield stringent constraints. However, 
it is suggestive that many of the objects of Nandra et 
al. (2007) "dream team" of relativistic AGN are NLSyls 
(Mkn 766, NGC 4051, NGC 5506, NGC 7314) or at least 
characterised by extreme variability on short (MCG-6-30-15, 
NGC 4395) or long (NGC 2992) timescales. Possible obser- 
vational biases are discussed (and ruled out) by Nandra et al. 
(2007). A trend towards relativistically broadened lines be- 
ing more common in Narrow Line AGN seems to be present 
in the stacked spectra of large samples of unobscured AGN 
( |Longinotti et al. 2008[ ). However, there doesn't seem to be 
a correlation between the detection of relativistically broad- 
ened iron lines and either the accretion rate or the black hole 
mass in the, e.g., FERO sample ( [Longinotti et al~2 008). 
Once again, the small size and incompleteness of the avail- 
able samples hamper firm conclusions. 
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